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The magnetic structures of small clusters of Fe, Mn, and Cr supported on a Cu(lll) surface have 
been studied with non-collinear first principles theory. Different geometries such as triangles, pyra- 
mids and wires are considered and the cluster sizes have been varied between two to ten atoms. The 
calculations have been performed using a real space linear muffin-tin orbital method (RS-LMTO- 
ASA). The Fe clusters are found to order ferromagnetically regardless of the cluster geometry. For 
Mn and Cr clusters, antiferromagnetic exchange interactions between nearest-neighbours are found 
to cause collinear antiferromagnetic ordering when the geometry allows it. If the antiferromagnetism 
is frustrated by the cluster geometry non-collinear ordering is found. A comparison between the 
calculated structures and ground states obtained from simplified Heisenberg Hamiltonians show that 
the exchange interaction varies for different atoms in the clusters as a result of the different local 
structure. 

PACS numbers: 75.75.+a,73.22.-f,75.10.-b 



I. INTRODUCTION 

The remarkable progress of experimental methodolo- 
gies with atomic resolution, such as the scanning elec- 
tron microscopy 1 (STM), has paved the way for stud- 
ies of nanoscale magnetic materials such as ad-atoms, 
clusters and wires deposited on surfaces. As a result of 
the reduced dimensions and symmetries for such systems, 
magnetic behaviour that differs from bulk materials can 
be found 2 -! 3 -. This attracts interest not only for the novel 
physics that can occur in these systems but also for the 
possibility to tailor the electronic and magnetic proper- 
ties by changing the structure and the local environment 
of the systems. 

Studies of systems consisting of only a few atoms can 
give valuable information on how the magnetic struc- 
ture evolves from single atoms towards the bulk be- 
haviour. Fe, Mn, and Cr are all known to exhibit in- 
teresting magnetic behaviour. While being ferromag- 
netic in the bcc phase, Fe in the fee phase has been 
found to exhibit a spin-spiral structure when synthe- 
sized as precipitates in a Cu matrix and calculations 
show that the magnetic structure is strongly dependent 
on the lattice parameter—. Cr has in bulk an incommen- 
surate antiferromagnetic spin density wave 6 - which can 
be tuned by creating superlattices with ferromagnetic or 
paramagnetic layers and varying the interface roughness 
and layer thickness 7 . When deposited on stepped sur- 
faces, Cr can be found to have non-collinear ordering 8 . 
Bulk Mn exhibits perhaps the most intriguing magnetic 
structure of all elements with a unit cell containing 58 
atoms 9 with a complex non-collinear antiferromagnetic 
magnetic structured Recently, several Mn based com- 
pounds, where the magnetic ordering is non-collinear due 
to geometric frustration between the magnetic moments 
of Mn atoms, have been studie d 11 ' 12 experimentally as 



well as theoretically. 

Free clusters of Fe, Mn, and Cr have been studied 
both experimentally and theoretically. Stern-Gerlach 
measurements on Fe^ 3 - clusters show ferromagnetic be- 
haviour while Mn clusters^ and Cr clusters^ show vary- 
ing small net deflections, which can be interpreted as 
the result of antiferromagnetic or even non-collinear mag- 
netic configurations. Calculations have shown that small 
Mn and Cr clusters can exhibit non-collinear magnetic 
orderin g 16 ' 17 : 18 in agreement with the Stern-Gerlach ex- 
periments. 

Supported transition metal clusters have also been 
studied extensively where a majority of the studies have 
been theoretical. Among the experimental studies, most 
work has been done on Fe, where Fe clusters deposited on 
a Ni surface have been found to be ferromagnetic with 
oscillating magnitude of the orbital moments— and Fe 
clusters supported on a graphite surface have been found 
to exhibit enhanced spin and orbital moments compared 
to bulk2°-. Among the theoretical studies, the reported 
calculations have mostly only considered collinear mag- 
netization densities. Monoatomic wires of Fe on Cu(lll) 
and Cu(001) show ferromagnetic behaviour with a strong 
magnetic anisotrop y 21 ' 22 . Small Mn clusters on Ag(001) 
have been found to exhibit magnetic bi-stabilit y 23 ' 24 
which is also the case for mixed clusters of FeMn and 
FeCr that have been found to have both ferro- and anti- 
ferromagnetic solutions close in energy- 2 ^. Early model 
calculations of supported equilateral triangular transi- 
tion metal clusters have shown that non-collinear order- 
ing can be obtained from the frustration due to antifer- 
romagnetic interactions between the cluster atomsj 26 ' 27 
Recent studies have found that small clusters of Mn and 
Cr become non-collinear when deposited on Ni(OOl) 28 
and Fe(OOl) 29 ' 30 surfaces due to competing exchange in- 
teractions between the cluster atoms and the surface 2 ^. 
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In a previous paper— we reported on non-collinear 
magnetic ordering for a selection of small Mn clusters 
supported on a Cu surface. In this paper we expand 
these results and present theoretical results concerning 
the magnetic ordering and interactions for Fe, Mn, and 
Cr clusters deposited on a Cu(lll) surface. The calcu- 
lations have been performed using the RS-LMTO-ASA 
method that is a first principles order-N method which 
has recently been extended to the treatment of non- 
collincar magnetism^. 

II. METHOD 

The RS-LMTO-ASA method is based on the LMTO- 
ASA technique 3 -^ and the Haydock recursion method 33 . 
The LMTO-ASA formalism provides an efficient, 
parameter-free, basis set for treating close packed metal- 
lic systems and the recursion method gives the ability 
to treat problems where translational symmetry is ab- 
sent and does also convey order-N scaling with respect 
to the number of nonequivalent atoms in the system. 
The recursion method does not directly solve the eigen- 
value problem as formulated in the DFT, but allows one 
to calculate the local density of states (LDOS) for the 
orbitals of the atoms in the selected system. The RS- 
LMTO-ASA method has successfully been used for a 
wide range of problems including bulk systems, multi- 
layers, embedded impurities and clusters and clusters on 
surfaces. Earlier and more detailed descriptions of the 
collinear implementation of the RS-LMTO-ASA can be 
found elsewher e - 34 ! 35 . 

In the local spin density approximation (LSDA)2&, 
the electron density is expressed through a 2x2 density 
matrix p which can be expressed in terms of the non- 
magnetic charge density n and the magnetization den- 
sity m as p = (nT + m ■ er)/2 where T is the 2x2 identity 
matrix and er = {cr x , a y , <r z } are the Pauli matrices. Self- 
consistent method o 37 ' 38 i 39 for calculating the electronic 
structure for non-collinear magnetization densities have 
existed for quite some time^P, and here we will focus on 
the specific details for treating non-collinear magnetiza- 
tion densities within the RS-LMTO-ASA. 

With the recursion method, the local density of states 
N(e) where e is the energy, is obtained as N(e) = 
— ^3trt/(e). Here Q(e) is the local Green's function 
Q(e) = (e — ri)~ l , where TC is the Hamiltonian. Simi- 
lar to the LDOS, the collinear magnetic density of states 
m(e) can be calculated as m(e) = — i3?tr(<7 z C?(e)). Since 
the Pauli spin matrix a z is diagonal in spin-space, the 
collinear magnetic density of states can be calculated us- 
ing only diagonal elements of the Green's function. If a 
generalized non-collinear magnetization density 

ro(e) = --3tr(ff0(e)) (1) 

where er = {<r x , cr y , a z } is sought, evaluation of the 
off-diagonal parts of the Green's function is in princi- 



ple needed. The off-diagonal elements of the Green's 
function are possible to obtain by performing the recur- 
sion starting from carefully selected linear combinations 
of muffin-tin orbitals 41 - or perform a computationally 
more demanding block recursion calculation 4 ^. However, 
in our implementation we avoid the evaluation of off- 
diagonal elements by applying successive unitary trans- 
formations U on the Hamiltonian, TC = UTLhfi . When 
the Hamiltonian is transformed in this way, the Green's 
function transform similarly; Q' — lAQU^ . 

Using the unitary property U^U = 1 and the fact that 
cyclic permutations of matrix multiplications conserve 
the trace of the product, the generalized magnetic den- 
sity of states m(e), can be written as 

m(e) = —$hr{oltfugiflU} = --3ir{er'£'}, (2) 

where er' is the Pauli matrices after the unitary trans- 
formation. The transformation matrix U is different for 
the three directions, and chosen so that, UajW = a' zl 
for j = x,y, z, to yield a diagonal representation. In the 
trivial case of j = z the unitary transformation is just 
the identity matrix. For the other directions, the unitary 
transformation corresponds to a spin rotation where U 
can be calculated using spin-i rotation matrices. De- 
composing the Hamiltonian into a spin-dependent part, 
B, and a spin-independent component, H , yields that U 
operates only on the spin-dependent part, 

TC = H + B -UctU ] . (3) 

From the transformed Hamiltonians, TC, the LDOS for 
the different directions can then be calculated using the 
regular recursion method and the magnetic density along 
the three directions can be obtained. From the three or- 
thogonal directions, the local magnetization axis is cal- 
culated and the LDOS for the local spin axis can be con- 
structed by taking the scalar product of the generalized 
magnetic density of states and the local magnetization 
vector. As all Hamiltonians are constructed within an ah 
initio LMTO-ASA formalism, all calculations are fully 
self-consistent, and the spin densities are treated within 
the local spin density approximation 3 ^. Since the recur- 
sion procedure is performed for three orthogonal direc- 
tions, the computational cost for each iteration is tripled 
compared with the collinear implementation of the RS- 
LMTO-ASA, but the linear scaling with respect to the 
number of nonequivalent atoms is retained. 

The calculations of the transition metal clusters have 
been performed by embedding the clusters as a pertur- 
bation on a self-consistently converged perfect Cu(lll) 
surface. The Cu surface has been calculated using the ex- 
perimental lattice parameter of Cu. As is usually the case 
for LMTO-ASA methods, the vacuum outside the surface 
needs to be simulated by having a number of layers of 
empty spheres above the Cu surface in order to provide 
a basis for the wave-function in the vacuum and to treat 
charge transfers correctly. After embedding the cluster 
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on the surface, the charge and magnetization densities 
of the cluster atoms and the neighboring Cu atoms and 
empty spheres are then recalculated until self-consistency 
is obtained while the electronic structure for atoms far 
from the cluster are kept unchanged to their unperturbed 
values. Structural relaxations have not been included in 
this study so the cluster sites have been placed on the 
regular fee lattice above the Cu surface. Earlier studied 
on supported transition metal clusters 4 ^ have shown that 
structural relaxations can change the magnetic properties 
of the clusters. On the other hand, in an experimental 
situation, small clusters as those considered in this study 
are usually constructed in an out-of-equilibrium situa- 
tion by manipulation with an STM tip and calculated 
equilibrium geometries might therefore not be relevant. 
The most relevant relaxation for these kind of artificially 
created clusters would be the distance between the clus- 
ter atoms and the substrate atoms and since a noble 
metal substrate is used in this study, the interaction be- 
tween clusters and substrate play a lesser role compared 
to the interactions between the cluster atoms. The clean 
Cu(lll) surface has been modeled by a large (>5000) 
slab of atoms and the continued fraction, that occurs 
in the recursion method, have been terminated with the 
Beer-Pettifor— terminator after 30 recursion levels. 

The non-collinear calculations have been performed 
without including the spin-orbit coupling. Since this 
term is neglected, a preferred spin axis does not exist 
in the system and the magnetic structures are thus only 
converged with respect to the directions of the magnetic 
moments relative to the other spin moments in the clus- 
ter. In order to minimize the risk of finding magnetic 
orderings that correspond to only a local minimum, sev- 
eral starting guesses where used for each system. 

The calculated magnetic structures can be analyzed 
in terms of the exchange interactions Jij between spins 
on atoms situated at sites i and site j. A well known 
connection between the exchange interactions and the 
magnetic ordering is given by the classical Heisenberg 
Hamiltonian, 
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Jijcos9ij, 
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where 9ij is the angle between the magnetic moment on 
site i and site j. Note that in Eqn. 01 the magnitude 
of the spins have been incorporated into the effective Jij 
interactions. 

In this work, we have calculated exchange interac- 
tions directly using the Liechtenstein formula 4 ^ as im- 
plemented in the RS-LMTO-ASA 4 ^ for a large selection 
of the considered clusters. The Jy's shown in this study, 
have been obtained from the ferromagnetic configuration 
of the clusters. Other magnetic configurations typically 
result in different values of the Jij 's, although the signs 
are seldom changed.— If the exchange interactions would 
be independent of 9ij, the magnetic structure could in 
principle be calculated by minimizing the Heisenberg 
Hamiltonian with J^'s calculated from a ferromagnetic 



configuration. This is not the case for the systems con- 
sidered in this work, which motivates a full non-collinear 
calculation of the magnetic structures. However, on a 
qualitative level the cause of the magnetic ordering, e.g. 
the effect of frustration or the competition between near- 
est and next-nearest interactions, can still be discussed 
in terms of the calculated exchange interactions. 

For a selection of Cr clusters, which are discussed 
in Sec. IIII CI our calculated magnetic structures has 
been compared to the structure found by minimizing 
the Heisenberg Hamiltonian for a fixed configuration of 
Jy-'s, where only nearest-neighbour interactions are fi- 
nite. The minimization of the Heisenberg Hamiltonian 
for these clusters has been performed by a genetic search 
algorithm 4 ^. 



III. RESULTS 
A. Fe clusters 

In Fig. Q] the magnetic structure of several Fe clusters 
are shown. Regardless of the geometry of the studied 
Fe clusters, we find the magnetic ordering in the clus- 
ters to always be ferromagnetic. The collinear magnetic 
structure for these clusters, which can be put in contrast 
with the non-collinear ordering found for fee structured 
Fe clusters embedded in bulk Cu^, could be caused by 
the fact that the decreased coordination of the surface 
clusters lead to a high-spin state which favours ferromag- 
netic coupling between neighbouring Fe atoms 49 . This is 
also consistent with an analysis by Lizarraga et al£Q. It 
should be pointed out here that large magnetic moments 
do not automatically lead to collinear magnetism. As 
we will see in the section below, Mn is an example where 
large moments result in an antifcrromagnetic interatomic 
exchange coupling, which on a frustrated geometry lead 
to non-collinear magnetism. In the case of Fe the large 
calculations result in large moments and a ferromagnetic 
interatomic exchange coupling. 

The spin moments for the Fe atoms in the clusters 
shown in Fig. [1] range between 3.45 hb for the atoms in 
the dimer to 2.56 hb for the central atom in the seven 

It has been shown 



atom cluster displayed in Fig. 1(d) 
for Co clusters in Cu(OOl) 35 and Fe clusters on Ni and 
Cu surfaces^ 1 -, that the magnetic moment has a linear 
behavior as function of the number of cluster neighbors 
around the site. The spin moments of the Fe clusters 
on Cu(lll) of Fig. [1] show a similar trend and depend 
almost linearly on the number of nearest Fe neighbours. 
For these ferromagnetic Fe clusters, the orbital moments 
were calculated and they where also found to depend on 
the number of nearest neighbours. The largest orbital 
moment was found to be 0.15 [Lb per atom for the atoms 
in the dimer, and the smallest orbital moment is 0.06 [Ib 
for the central atom in the cluster shown in Fig. l(d)| 

Due to the strong correlation between the magnetic 
moment for the atoms in the Fe clusters and the num- 
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FIG. 1: The geometries for clusters of Fe atoms on a Cu(lll) 
surface. All Fe clusters are found to exhibit a ferromagnetic 
ground state regardless of the cluster structure. The num- 
bers indicate the atom projected spin moment of the different 
atoms. 
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ber of Fe neighbours, the magnetic moments obtained 
above can in principle be used to predict the total mag- 
netic moment of any Fe cluster as long as the shape is 
determined and the cluster is planar. This would indi- 
cate that the magnetic moment per atom for a perfect 
monolayer of Fe atoms on a Cu(lll) surface would be 
~ 2.7 [is which is in good agreement with earlier calcu- 
lations of Fe monolayers on C u 52 ' 53 i 54 ' 55 . Our findings of 
ferromagnetic coupling indicate that a single Fe mono- 
layer would be ferromagnetic in agreement with Ref. 
whereas Ref. [HH found that a single row antiferromagnetic 
order would be the most stable magnetic configuration. 
This discrepancy may be explained by the use of different 
lattice parameters in Refs. [52| and [53|. 



B. Mn clusters 

In a previous paper~ we showed that due to antifer- 
romagnetic coupling between nearest neighbour atoms in 
Mn clusters deposited on Cu one finds either a collincar 
antiferromagnetic structure or, if frustration occurs due 
to the cluster geometry, a non-collinear magnetic struc- 
ture. A collection of frustrated cluster geometries with 
triangular shapes is shown in Fig. [2] The magnetic mo- 
ments obtained for each atom of the studied Mn cluster- 
are shown in Table. [T] For the equilateral triangle in 



Fig. 2(a) a non-collinear arrangement with angle of 120° 



FIG. 2: The calculated magnetic ordering for triangular Mn 
clusters on a Cu(lll) surface. For all geometries except the 
isosceles triangle shown in (b) and the six atom triangle in (c), 
non-collinear solutions are obtained due to magnetic frustra- 
tion. 



non-collinear solution and a frustrated collinear antifer- 
romagnetic solution, with two spins parallel to each other 
and anti-parallel to the third spin is 13 meV/atom. The 
ferromagnetic solution was found to have an energy of 
102 meV/atom higher than the stable non- collin ear so- 
lution. The isosceles triangle shown in Fig. |2(b)| has an 
antiferromagnetic collinear ground state which indicates 
that the exchange coupling between the two Mn atoms 
furthest from each other is either very small compared to 
the antiferromagnetic nearest-neighbour exchange cou- 
pling or has the opposite sign (i.e. ferromagnetic). As 
the s ize of the triangular clusters increases (Figs. 2(c) to 
2(e) I, the behaviour becomes more intricate. 



Although 

the six-atom triangle (Fig. 2(c)| by the analogous ge- 
ometry as the three atom triangle in Fig. 2(a) could be 



between the magnetic moments is the most stable config- 
uration. The total energy difference between the stable 



expected to align in a structure with 120° between neigh- 
bouring atoms, it is in fact a collinear antiferromagnetic 
order that is the most stable solution. The energy differ- 
ence between the collinear structure and a non-collinear 
structure was found to be 18 meV/atom. The cause for 
the preferred collinear order is the different environment 
for the corner atoms who only have two nearest neigh- 
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bours compared to the three central atoms who have four 
nearest neighbours each. The reduced coordination for 
the corner atoms causes their antiferromagnetic exchange 
coupling to nearest neighbours to be enhanced. The cal- 
culated JijS confirm this behaviour since the strength of 
the exchange interaction between a corner atom and a 
nearest-neighbour amounts to -27 meV while it is -12 
meV between two central atoms. A similar mechanism 
can be expected for the nine-atom cluster displayed in 



Fig. 2(d) but for this geometry the atoms that are not 
situated at the corners of the triangle have three nearest 
neighbours due to the hole in the middle of the cluster. 
Therefore the difference in the local geometry is smaller 
between the corner atoms and the central atoms which 
leads to a more delicate balance of the exchange cou- 
plings. The resulting structure has the moments point- 
ing in three different directions instead of two directions 
which would be the case for a collinear antiferromagnetic 
solution. The angle between two neighbouring central 
atoms is 152° while the angle between a corner atom and 
a nearest- neighbour is 104°. Contrary to what was found 
for the six atom cluster in Fig. 2(c) it thus appears that 



the exchange coupling is larger between central atoms 
then between a corner atom and a central atom. This be- 
haviour is supported by the calculated Jys where the cou- 
pling between an corner atom and a nearest-neighbour is 
-9meV while the exchange coupling is found to be -36 
meV between two central atoms. The non-collinear solu- 
tion for the nine atom cluster has a total energy which 
is 13me V/atom lower than a collinear antiferromagnetic 
solution. 

The analysis of the final ten-atom triangle shown in 
Fig. 2(e) is even more complicated. From a geometrical 
view, this cluster has three nonequivalent sites; the three 
corner atoms, the six atoms neighbouring to the corner 
atoms and the central atom. The magnetic structure does 
however have a lower symmetry that can be described by 
decomposing the cluster into the central atom and three 
'sub-triangles', consisting of the three atoms closest to 
each corner of the cluster. Within each sub-triangle, the 
three atoms couple to each other in a geometry that re- 
sembles the 120° structure of a single three atom trian- 
gular cluster, but since exchange interactions from other 
neighbouring atoms are present as well, the angles be- 
tween the moments in these sub-triangles vary between 
146° and 104°. All angles between the moments for the 
atoms in the ten atom cluster can be seen in Table|TTJ The 
influence of the different number of neighbours for the 
cluster atoms determine their magnetic moments where 
the corner atoms have a magnetic moment of 4.3 /ib, the 
central atom has 2.6 /j,b and the magnetic moment for 
the remaining six atoms is 3.6 /zb- 

Atomic wires constitute a group of nanostructures that 
has attracted a lot of attentio n 3 i 56 i 57 . We have calcu- 
lated the magnetic structure for wires of Mn atoms with 
lengths between two to nine atoms. The wires are ori- 
ented along a 110 direction on the Cu surface. The total 
energy differences, per cluster atom, between the anti- 



TABLE I: Magnetic moments (in ps) for atoms in the clus- 
ters displayed in Fig. 2(a)|2(e) The atoms are numbered 
in the left column, starting from the leftmost atom, increas- 
ing around the cluster in the clockwise direction and, for the 
largest cluster, ending with the central atom. 
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gles between moments for atoms in the cluster displayed in 



Fig. 2(e) The atoms are numbered starting from the leftmost 
atom, increasing around the cluster in the clockwise direction 
and ending with the central atom. 



Atom 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 





146 


52 


147 


57 


20 


90 


150 


104 


103 


2 


146 





98 


52 


155 


158 


56 


60 


107 


52 


3 


52 


98 





146 


107 


60 


56 


158 


155 


52 


4 


147 


52 


146 





104 


150 


90 


20 


57 


103 


5 


57 


155 


107 


104 





47 


138 


96 


49 


150 


6 


20 


158 


60 


150 


47 





108 


141 


96 


107 


7 


90 


56 


56 


90 


138 


108 





108 


138 


53 


8 


150 


60 


158 


20 


96 


141 


108 





47 


107 


9 


104 


107 


155 


57 


49 


96 


138 


47 





150 


10 


103 


52 


52 


103 


150 


107 


53 


107 


150 






ferromagnetic and the ferromagnetic magnetic configu- 
rations for the Mn wires are shown in Fig. [3] A large 
energy difference of 96 meV/atom is found for the dimer 
whereas the energy differences for the longer wires are 
significantly smaller. If only nearest neighbour interac- 
tions played role one would expect an energy difference 
(per atom) between ferromagnetic and antiferromagnetic 
coupling with a functional form J(l--i-), where N is the 
number of atoms in the wire. Hence for long chains this 
energy difference should be equal to J whereas the energy 
difference would continuously become smaller and reach 
the value J/2 for N=2. The data in Fig. [3] does not dis- 
play this trend, which suggests that next-nearest neigh- 
bour interactions are important and/or that the value of 
J depends on the number of atoms of the cluster. 

For the dimer and trimer only a collinear antiferro- 
magnetic solution is found, whereas, for longer chains a 
slightly canted non-collinear order is also found. For al- 
most all of the longer chains, the non-collinear solutions 
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TABLE III: Calculated exchange parameters Jy (in meV) for 
the Mn pentamer shown in Fig. 3] The atoms are numbered 
from left to right. 
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FIG. 3: Total energy differences between antiferromagnetic, 
E AF , and ferromagnetic, E FM , configurations of wires of Mn 
atoms, oriented along a (110) direction on a Cu(lll) surface. 




FIG. 4: Calculated magnetic configurations for a five atomic 
Mn wire, oriented along a (110) direction on a Cu(lll) sur- 
face. 



are unstable, but resembles the collinear antiferromag- 
netic solution closely, both in energy, in all cases less than 
0.5 meV/atom, and in the angular difference, where the 
deviation from the collinear structure is smaller than 3° 
per atomic pair. 

The Mn pentamer is however an interesting exception 
from the behavior of the other wires, and for this system, 
the ground state is actually found to be a non-collinear 
configuration, shown in Fig. 3] The non-collinear config- 
uration for the pentamer can be described by the angle 
between an edge atom and its nearest neighbour, which 
is 170°, and the angle between two neighbouring central 
atoms, which is 155°. The energy difference between the 
non-collinear and the antiferromagnetic solutions for the 
pentamer is 2 meV/atom. 

The straight wires does not have geometries that cause 
frustration in the same way as the triangular clusters 
mentioned earlier, so a probable cause for the non- 
collincarity of the Mn pentamer is the competition be- 
tween ferromagnetic and antiferromagnetic exchange in- 
teractions between the different atoms in the cluster. 
Since the nearest-neighbour interactions are always an- 
tiferromagnetic for the Mn clusters (at least for the 
nearest- neighbour distance used in this study), more 



long-range interactions must play a role in destabiliz- 
ing the collinear magnetic state. In order to examine 
the size and range of the exchange interactions, we have 
calculated exchange coupling parameters J+j for the five 
atom wire shown in Fig. 3J The values of the exchange 
parameters are shown in Table IIIII where the i and j are 
chosen so that site 1 and 5 are the edge atoms and site 
3 is the central atom. Hence site 2 is nearest neighbour 
to site 1 and 3 and the nearest neighbours for site 4 is 
site 3 and 5. It may be observed that the exchange in- 
teractions are strongest and antiferromagnetic between 
nearest neighbours, with a smaller long range interaction 
that oscillates between ferromagnetic and antiferromag- 
netic coupling. The largest magnitude for the exchange 
interaction is obtained for J12 and J45, i.e. between an 
edge atom and its nearest neighbour. Furthermore, Ta- 
ble [HT] shows that although the nearest neighbour inter- 
actions have the largest magnitude, the more long ranged 
interactions always seem to counteract the nearest neigh- 
bour interactions. This might not be obvious from the 
values in Table HTT1 but as a clarifying example we can ex- 
amine the exchange interactions between atom 1 and the 
other atoms. The negative nearest neighbour interactions 
in the pentamer would prefer an antiferromagnetic order 
so that atom 1 would be ferromagnetically coupled to 
atom 3 and atom 5, and antiferromagnetically coupled to 
atoms 2 and 4. However, the calculated exchange inter- 
actions in Table IIIII show that the J13 and J15 are in fact 
negative while J14 is positive, thus competing against the 
antiferromagnetic nearest-neighbour ordering, which re- 
sults in the non-collinear magnetic state shown in Fig. 2] 



C. Cr clusters 

In Fig. [5] the geometries and calculated magnetic con- 
figurations for a selection of Cr clusters are shown. Cal- 
culations for Mn clusters with similar geometries can be 
found elsewhere^ The magnetic moments obtained for 
each atom of the Cr clusters shown in Fig.[5]can be seen in 
Table. IIVI The magnetic structures are for most of the 
clusters quite similar to the calculated magnetic struc- 
ture for Mn clusters although certain differences occur, 
as will be commented on below. Both the Cr dimer and a 
straight trimer (not shown) orders antiferromagnetically 
while a three atom triangle has the same non-collinear 
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(a) 



(b) 




(e) 



(f) 



FIG. 5: The calculated magnetic ground state for Cr clusters 
on a Cu(lll) surface. 



TABLE IV: Magnetic moments (in [ib) for atoms in the clus- 
ters displayed in Fig. 5(a)|5(f) The atoms are numbered 
starting from the leftmost atom, counting around the cluster 
in the clockwise direction and, if applicable, ending with the 
central atom. 





5(a) 


5(b) 


5(c) 


5(d) 


5(e) 


5(f) 


1 


3.99 


3.99 


3.82 


3.54 


3.87 


3.57 


2 


3.99 


3.89 


3.82 


3.54 


3.65 


3.28 


3 


3.99 


3.92 


3.82 


3.54 


3.89 


3.89 


4 


3.50 


3.89 


3.82 


3.54 


3.65 


3.28 


5 




3.99 


3.82 


3.54 


3.87 


3.57 


6 






3.82 


3.54 


3.65 


3.28 


7 








3.12 


3.89 


3.89 


8 










3.65 


3.28 


9 












2.81 



structure as its Mn counterpart, as shown in Fig. 2(a) 
It can be noted that this non-collinear structure has also 
been reported from calculations on Cr clusters, with the 
same geometry, supported on Au(lll) i ' The cluster in 
Fig. 5(a) has a collinear antiferromagnetic ground state 
since the edge atoms only have one nearest neighbour 
and therefore no frustration occurs. The pentamer in 
Fig. |5(b)j also exhibits an antiferromagnetic ground-state, 
which in contrast to the non-collinear behaviour of the 
Mn pentamer, is purely collinear. This finding indicates 



TABLE V: Calculated exchange parameters Jij (in meV) for 
the Cr pentamer shown in Fig. 5(b) The atoms are numbered 
from left to right. 



i\j 


1 


2 


3 


4 


5 


1 




-143 


4.6 


18.9 


15.3 


2 


-143 




-97.0 


-40.5 


18.9 


3 


4.6 


-97.0 




-97.0 


4.6 


4 


18.9 


-40.5 


-97.0 




-143 


5 


15.3 


18.9 


4.6 


-143 





that the magnetic structures of the Cr clusters are more 
strongly dependant on the nearest-neighbour exchange 
coupling than the Mn clusters are. The calculated ex- 
change parameters Jij for the Cr pentamer is shown in 
Table [V] Compared with the exchange interactions for 
the Mn pentamer in Table IIIII we see that the nearest- 
neighbour interactions are indeed larger between the Cr 
atoms. On the other hand, the more long ranged inter- 
actions also have larger magnitudes in the Cr pentamer 
than for the Mn counterpart. The exchange interactions 
between atoms further from each other do however not al- 
ways compete against the nearest-neighbour interactions 
as was the case for the Mn pentamer. 

The collinear antiferromagnetic behaviour found for 
the pentamer also occurs for the six atom ring displayed 
in Fig. 5(c)[ which is expected since the geometry does 
not cause frustration for the nearest-neighbour interac- 
tions. The cluster in Fig. 5(d) has a symmetric non- 
collinear ground state with an angle between two neigh- 
bouring atoms on the rim of the cluster of 157° and be- 
tween the central atom and any outer atom the angle 
is 101°. This can be compared with the energy min- 
imum obtained when minimizing the nearest-neighbour 



Heisenberg Hamiltonian where all nearest neighbour Jij 's 
are set to be negative but equal. The Heisenberg model 
would give equilibrium angles of the same cluster ge- 
ometry of between 151° outer neighbours and 104° be- 
tween the central atom and any neighbour. Although the 
agreement between our calculated ground state and the 
Heisenberg minimum is good, it should be noted that 
due to the difference in the local structure around the 
central and outer atoms the exchange parameters Jij 
should be different between two outer atoms compared 
to J^ 's connecting to the central atom. This difference 
can be considered in a simple model analysis by damp- 
ing the strength of the exchange parameters where the 
central atom is connected, and for a damping of 20% 
for these exchange parameters, the Heisenberg Hamilto- 
nian approach yields an energy minimum with the angles 
of 157° for neighbouring outer atoms and 101° between 
the central atom and an edge atom which are in perfect 
agreement with our calculated angles. 



Since the atoms of the hollow cluster in Fig. 5(e) do not 
all have only two nearest neighbours, as is the case for 
the other ring like geometry of Fig. 5(c) the cluster can 



not have an unfrustrated antiferromagnetic solution. De- 
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TABLE VI: Angles between magnetic moments for atoms 
in the cluster displayed in Fig. 5(e) The atoms are num- 
bered starting from the leftmost atom and increasing along 
the clockwise direction of the cluster. 



;om 


1 


2 


3 


4 


5 


6 


7 


8 


1 





169 


56 


56 


112 


56 


56 


169 


2 


169 





112 


133 


56 


133 


112 





3 


56 


112 





112 


56 


112 





112 


4 


56 


133 


112 





169 





112 


133 


5 


112 


56 


56 


169 





169 


56 


56 


6 


56 


133 


112 





169 





112 


133 


7 


56 


112 





112 


56 


112 





112 


8 


169 





112 


133 


56 


133 


112 







(a) 



(b) 



(c) 



FIG. 6: The calculated magnetic ordering for pyramid shaped 
clusters on a Cu(lll) surface. Fig. 6(a) shows a Fe cluster 
with a ferromagnetic solution. Fig. 6(b) shows a Mn cluster 
and Fig. 6(c) shows the Cr pyramid. 



TABLE VII: Angles between magnetic moments for atoms in 
the cluster displayed in Fig. 5(f) The atoms are numbered 



starting from the leftmost atom, increasing along the clock- 
wise direction of the cluster and ending with the central atom. 



Atom 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 





156 


40 


95 


74 


95 


40 


156 


69 


2 


156 





137 


86 


95 


86 


137 





128 


3 


40 


137 





137 


40 


137 





137 


40 


4 


95 


86 


137 





156 





137 


86 


128 


5 


74 


95 


40 


156 





156 


40 


95 


69 


6 


95 


86 


137 





156 





137 


86 


128 


7 


40 


137 





137 


40 


137 





137 


40 


8 


156 





137 


86 


95 


86 


137 





128 


9 


69 


128 


40 


128 


69 


128 


40 


128 






scribing the cluster with a Heisenberg Hamiltonian with 
equal and negative Jy 's would yield a ground state with 

As seen in Fig, 



120° between neighbouring atoms. As seen in Fig. |5(e) 
our calculated magnetic structure differs from the Heisen- 
berg minimum, it is instead described with three differ- 
ent nearest-neighbour angles. The upper and lower edge 
atoms have an angle of 112° to their neighbours while 
the leftmost and rightmost of the atoms have an angle 
of 169° to their nearest neighbours. The third angle is 
that between two atoms with three neighbours each and 
they have an angle of 133° between them. This struc- 
ture can be explained in a similar way as the cluster in 
Fig. 5(d) could, with different local geometries resulting 
in different strengths of the exchange coupling and thus 
different J^'s for different atoms. The angles between 
the different magnetic moments for the cluster shown in 
Fig. 5(e) are given in Table. I VII 



Filling the empty site in the middle of the cluster in 
Fig. 5(e)| gives the cluster geometry shown in Fig. 5(f) 



This additional atom causes the magnetic structure to 
be even more complex. Since the symmetry is lowered 
compared to the cluster in Fig. 5(d) the central atom 
does not have the same angle towards all of its neigh- 
bours. Instead two angles are needed to describe the 



5(f) are given in 



structure of the neighbours of the central atom. One of 
these is the angle of 69° which the central atom makes 
towards the leftmost and rightmost atoms and the other 
angle connects the central atom with the remaining four 
neighbours and the size of this angle is 128°. The up- 
per and lower edge atoms makes and angle of 137° with 
their neighbours. The angles between different magnetic 
moments for the cluster shown in Fig. 
Table. IVIll 



D. Three dimensional clusters 

So far the studied clusters have all been confined in 
one layer above the Cu surface. However, our method 
can treat three dimensional clusters as well. In order to 
demonstrate this we show the obtained magnetic config- 
urations for a pyramid-like tetrahedron shaped cluster 
in Fig. As expected from the results for Fe clusters 
reported earlier in this work, the Fe cluster (Fig. 6(a) ) 



exhibits a ferromagnetic order. The atom situated on top 
of the pyramid has a magnetic moment of 3.40/iBwhile 
the three Fe atoms closer to the Cu surface have a mag- 
netic moment of 3.11/xb- For the Mn cluster shown in 
Fig. |6(b)| a non-collinc ar stru cture is found and for the Cr 
pyramid, shown in Fig. 6(c) a collinear antiferromagnetic 
solution is found. A model Heisenberg Hamiltonian, as 
in Eqn. [4j with only antiferromagnetic nearest-neighbour 
exchange parameters, Jy , of equal size, yields a two-fold 
degenerate ground state, either a collinear antiferromag- 
net or a non-collinear tetragonal configuration with 109° 
between neighbouring angles. The calculation for the Mn 
pyramid over Cu(lll) shows angles which are slightly dis- 
torted relative to those of the free pyramid, around 116 
degrees between the base atoms close to the substrate 
and angles of about 100 degrees between the base site 
and the top one. It is peculiar that the Mn cluster has 
the non-collinear tetragonal configuration as the ground 
state whereas the Cr cluster has the collinear antiferro- 
magnetic ground state solution. A possible explanation 
for the non-collinear ground state of the Mn pyramid 
could be that, similar to the situation for several of the 
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planar clusters, the reduced neighbour coordination for 
the top atom compared to the atom in the base of the 
pyramid yields different exchange interaction strengths 
between the atoms in the cluster. The calculated ex- 
change parameters confirm this since the exchange inter- 
action between a top and a base atom is —83 meV while 
the interaction between two base atoms —46 meV. The 
reduced neighbour coordination also affects the magnetic 
moment for the top atom which is 4.5 /is compared to 
4.0 fiB for the base atoms. However, the situation with 
different exchange parameters is also present for the Cr 
pyramid where the interaction strength between the top 
atom and a base atom is —171 meV compared to —98 
meV for the coupling between two base atoms. The mag- 
netic moment for the Cr atom on top of the pyramid is 
4.2 \xb while the magnetic moment for the base atoms is 
3.6 hb- The difference in the magnetic ordering found 
for the Cr and Mn pyramids indicates that the bilinear 
exchange terms J cannot always describe the magnetic 
interactions between the atoms in supported magnetic 
clusters, a fact which previously been suggested for mag- 
netic dimers on surfaces.— It can be noted that the differ- 
ence in the total energy between the non-collinear ground 
state structure and the antiferromagnetic solution for the 
Mn cluster is 25 meV per atom while the corresponding 
difference for the Cr pyramid is -15 meV per atom. 

IV. CONCLUSIONS 

We have studied the magnetic structure of small clus- 
ters of Fe, Mn, and Cr supported on a Cu(lll) sur- 
face with non-collinear, first principles calculations. The 
studied Fe clusters are found to order ferromagnetically 



regardless of the cluster geometry. For Mn and Cr clus- 
ters, antiferromagnetic exchange interactions between 
nearest-neighbours are found to cause either collinear an- 
tiferromagnetic ordering or non-collinear ordering. The 
non-collinear ordering occurs when the cluster geometry 
is such that an antiferromagnetic arrangement becomes 
frustrated. The calculations have been accompanied by 
comparisons with calculated effective exchange interac- 
tions as well as with ground states obtained from a simpli- 
fied Heisenberg Hamiltonian and the comparisons show 
that the exchange interactions vary for different atoms 
in the clusters as a result of the different local struc- 
ture. Differences between the magnetic ordering for Mn 
and Cr clusters are found where Cr clusters seem to pre- 
fer collinear solutions to a higher degree while Mn clus- 
ters can exhibit non-collinear configurations even for un- 
frustrated cluster geometries. Comparisons with model 
Hamiltonians show that the magnetic structure of certain 
clusters can be explained by a simple nearest-neighbour 
Heisenberg Hamiltonian while other cluster geometries 
cause more complex behaviours. 
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